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Abstract. With a simple power-law approximation of high-redshift (> 3.5) star formation 
history, i.e., p*(z) oc [(1 + z)/4.5] _Q , we investigate the reionization of intergalactic medium 
(IGM) and the consequent Thomson scattering optical depth for cosmic microwave back- 
ground (CMB) photons. A constraint on the evolution index a is derived from the CMB 
optical depth measured by the Wilkinson Microwave Anisotropy Probe (WMAP) experiment, 
which reads a ~ 2.18 lg ^4C — 3.89, where the free parameter jK, is the number of the escaped 
ionizing ultraviolet photons per baryon. At the same time, the redshift zj at which the IGM 
is fully ionized can also be expressed as a function of a as well as jV^. By further taking 
into account the implication of the Gunn-Peterson trough observations to quasars for the 
full reionization redshift, i.e., 6 < zr < 7, we obtain 0.3 < a < 1.3 and 80 < < 230. 
For a typical number of ~ 4000 of ionizing photons released per baryon of normal stars, the 
fraction of these photons escaping from the stars, / esc , can be constrained to within the range 
of (2.0 - 5.8)%. 
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1 Introduction 

The cosmic history of star formation (especially at high redshifts) is one of long outstanding 
problems in astrophysical cosmology, which is being probed by a growing number of obser- 
vations such as through Lyman break galaxies (LBG) [1, 2], Ly-a emitter galaixes [3], and 
gamma-ray bursts (GRBs) [4-8]. One of the most important consequences of star forma- 
tion is the reionization of intergalactic medium (IGM) that has been recombined at redshift 
z oj 1100, although some other sources (e.g., accreting massive black holes) may also provide 
contributions to the reionization. The answer to the question that how and when IGM was 
reionized must hold many clues to the nature of the first generation of light sources and also 
to the onset of structure formation in cold dark matter cosmologies. During reionization, the 
librated free electrons can scatter cosmic microwave background (CMB) photons, resulting 
in an increase in the degree of the anisotropies of both CMB temperature and polarization 



Therefore, since the launch of the Wilkinson Microwave Anisotropy Probe ( WMAP) 
satellite, it has become of great interest to derive constraints on the cosmic histories of star 
formation and reionization from CMB observations. As one of the most direct results, the 
WMAP- 7 observation infers a Thomson scattering optical depth of CMB photons r = 0.088± 
0.015 [10, 11]. Then, the redshift for full reionization can easily be derived to zt = 10.5 ± 1.2 
with the assumption that the transition from a fully neutral to a fully ionized IGM is abrupt. 
However, deeper implications of the optical depth always invoke an elaborate description 
for the star formation history, either by theoretical modelings in the hierarchical formation 
model [12, 13] or by observational determinations from some specific observations [14]. In 
this paper, we alternatively take a power-law assumption for the high-redshift star formation 
history as an effective description, where all of the relevant physical details are packed into 
a free parameter (i.e., the power law index). The power-law behavior could be regarded as 
the first order approximation of the realistic history, which in principle allows us to get some 
general and substantial inferences. 

2 Model 

First of all, all analyses throughout the paper is based on the assumption that cosmic reion- 
ization is dominated by the early generation of stars. The rate of ionizing ultraviolet photons 
escaping from stars into IGM reads 



[9]. 
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where p*(z) is the star formation rate (SFR), is the mass of a baryon, iV 7 is the number 
of ionizing ultraviolet photons released per baryon of the stars, and / esc is the fraction of 
these photons escaping from the stars. In literature, the value of iV 7 is usually taken to be 
~ 4000 for a Salpeter stellar initial mass function and a metallicity ~ 0.05Z Q (e.g., [15]). 
For the value of / csc , a wild range from a few per cent to ~ 20 per cent was suggested by 
some observations [16-19]. Similar values were also predicted by some theoretical modelings 
and simulations. For example, Razoumov & Sommer-Larsen [20] found that f esc evolves 
from ~ 1 — 2 percent at z = 2.39 to ~ 6 — 10 percent at z = 3.6. In this paper, we would 
use the product of iV 7 and / C sc as a single parameter denoted by JY^ = iV 7 / esc , which is 
considered to be completely free. Accompanying with star formation, IGM was reionized. 
The density of the liberated electrons n e can be connected to the hydrogen density nn by 
n e = (1 + y)riHX, where x = nu n/nu represents the fraction of ionized hydrogen and y is 
introduced by considering the ionization of helium. By assuming that the helium was only 
once ionized, we have y = Y/{AX) ~ 0.08, where X = 0.7523 and Y = 0.2477 are the mass 
fractions of the hydrogen and helium, respectively, at the reionization era [21]. Following ref. 
[15], the evolution of x can be determined by 



dx 
dz 



(i + yK 



a B C(l + zf(l + y)n^x 



dz 



where = 1.9 x 10~ 7 cm" 3 is the present number density of hydrogen, 0.3 = 2.6 x 
10~ 13 cm'^s" 1 is the recombination coefficient of netral hydrogen for electron temperature 
of ~ 10 4 K, C = (n^ n )/(^H n ) 2 i s the clumping factor of ionized gas, and dt/dz = 
-l/[H(z)(l + z)\ with H(z) = H ov /n m (l + z) 3 + n A , H = 71 km s^Mpc^ 1 , n m = 0.267, 
and Ov = 0.734 [10]. 

For a given reionization history x(z), the CMB optical depth can be calculated by 
integrating the electron density times the Thomson cross section along proper length 

r = J OTn e dl 

= -(l + y)a T n^c (1 + zfx{z)— dz. (2.3) 
Jo &Z 

The WMAP- 7 observation gave r = 0.088 ± 0.015, which hence provides a robust constraint 
on the reionization history. There seems a trouble in the precise definition of the upper limit 
of the above integral, because of the present poor knowledge of the onset of star formation. 
However, due to the rapid decease of the degree of ionization with increasing redshifts (see the 
results presented in the next section), the CMB optical depth could be mainly contributed 
by the electrons at relatively low redshifts z z^. More importantly, since the WMAP 
experiment is probably insensitive to too high redshifts [10, 28], the value of Zh actually can 
not be taken to extremely high. Therefore, when we solve eq. (2.2) from high to low redshifts, 
we somewhat arbitrarily take the initial redshift z^ to be a moderate value of 20 at which the 
transition from Population (Pop) III to I /II stars 1 may take place [23-26]. Correspondingly, 
the value of x at z^ is in principle determined by both the residual electrons left from the 



1 The Pop III-I/II transition is usually considered to take place when the cosmic metallicity is enriched to 
~ 1O~ 3 ' 5 Z0. However, the recently discovered Caffau star showed that the critical metallicity could be much 
smaller due to dust cooling, i.e., ~ 10 -5 ' Zq [22]. So the transition could proceed more rapidly than previously 
considered. 
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Figure 1. Star formation history with 0.3 < a < 1.3 (shaded band), which is constrained by both 
the WMAP-7 CMB optical depth (i.e., r = 0.088) and the Gunn-Peterson trough observations (i.e., 
6 < Zttitt ^5 7). Some observationally-determined SFRs are also presented for comparisons. 

recombination era [27] and, in particular, the electrons ionized by possible Pop III stars at 
z > Zfi- For a generally consideration, a wide range of < x(zf l ) < 0.1 would be taken into 
account in our calculations. Here all of the uncertain physics above redshift Zh is packed into 
the free parameter x(zh)- Finally, it should be noticed that the contribution to the CMB 
optical depth by the residual electrons was already accounted for when deriving r from the 
CMB data. Therefore, in principle, the fraction of these electrons needs to be deducted from 
the total x(z) when we calculate the optical depth by eq. (2.3). Nevertheless, in view of the 
actually small fraction of the residual electrons on the order of 10 -3 [27], such a deduction 
would be ignored in our calculations. 

Solving eqs. (2.2) and (2.3) depends on two crucial but controversial astrophysical 
inputs, i.e., the star formation rate p*(z) and the clumping factor C(z). Following a series 
of measurements of SFRs, especially the complication of Hopkins & Beacom [29] which is 
shown in figure 1 (solid circles), a consensus on the history of star formation now emerges 
up to redshift z ~ 3.5. The behavior includes a steady increase of star formation from z = 
to z = 1, and a following plateau up to z ~ 3.5, which can be empirically expressed by 



with a local rate p*(0) = 0.02 M Q yr -1 Mpc -3 . For z > 3.5, a clue to a decrease in the 
SFRs was seemingly implied by the Hopkins & Beacom's data, which was further confirmed 
by the observations of LBGs and GRBs. Anyway, in view of the obvious insufficiency of the 
high-redshift observations, the tendency of the star formation history to high redshifts is in 
fact still ambiguous. Therefore, all of the cases including those in which the history continues 
to plateau, drops off, or even increases were taken into account in previous literature (e.g., 
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[30]). Instead of a theoretical modeling or observational determination, we here introduce a 
free parameter a to parameterize the high-redshift history as a power law 

p,(z) =0.2(1^) M yr^Gpc" 3 , for z > 3.5. (2.5) 

Such a simple power-law behavior would be regarded as the first order approximation of 
the realistic history. Finally, for the clumping factor, its value can be estimated by numer- 
ical simulations (e.g., [31, 32]) and by semianalytical models (e.g., [33-36]). All of these 
calculations showed that the clumping factor decreases with increasing redshifts, which can 
approximately be described by a power law. Following a recent simulation given by ref. [32], 
here we adopt 

C(z) = 2.9 0-^) > ^ z > 5 (2.6) 
which has a fixed C(z < 5) = 2.9. 



3 Results 

Substituting eqs. (2.4-2.6) into eq. (2.2), a cosmic reionization history, x(z), can be solved 
for a set of values of the parameters JV^ and a. By equating the CMB optical depth for this 
reionization to the WMAP-7 one (i.e., r = 0.088), we get a nearly linear relationship between 
the parameters Ig^y and a, as presented in the top panel of figure 2, which can be fitted by 

a = alg^-6 (3.1) 

with a = 2.18 ± 0.05 and b = 3.89 ± 0.01. Here the error bars of the coefficients are resulted 
from the uncertainty of the parameter x(zh)- Moreover, the small values of the error bars 
indicate that the jY^ — a relationship is actually insensitive to the value of x(zh), which 
demonstrates that it is acceptable to ignore the deduction of the fraction of residual electrons. 
With the above relationship, we present the solutions to eq. (2.2) in figure 3 for different 
values of a, where the redshift at which the IGM is fully ionized (denoted by zj) can be 
found. Subsequently, we can obtain the relationship between zj and JKj, as shown in the 
bottom panel of figure 2. 

As also shown by figure 3, the fittings to the reionization histories for Zf < z < 10 
exhibit a power-law behavior as 

/ 1 + z \ -( Q+1 - 7 ) 

x ra {TT7 f J ' for z ~ Zf - (3,2) 

Substituting the above approximation into eq. (2.3), an analytical expression for the optical 
depth can be obtained 

2a T n° H c(l + y) 

T ~ — 



(1 + ^/) 3/2 ( 1 + -^-) 
K V \ 2 + lOa/ 



n 1/2 



(3.3) 



where 77 = 1 — [(1 + Zh)/(1 + Zf)\~ a ~ 1 ^ and the parameter Q\ is omitted due to £l\ <C 
fJ m (l + zj) 3 . For the WMAP- 7 optical depth r = 0.088, the redshift for full reionization as 
a function of a can be expressed 

z ~ 11.7 ( 1 + ~ 2/3 _ 1 (3 4) 

t V 2 + 10a/ K ' 
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Figure 2. Relationships required by the WMAP-7 CMB optical depth (t = 0.088) between the 
parameters JY^ and a (top panel) or jY^ and Z{ u n (bottom panel), as shown by the solid bands. 
The width of the bands is determined by the adopted range of the initial degree of ionization as 
< x(zh) < 0.1. The possible range of Zf denoted by the hatched band is inferred from the Gunn- 
Peterson trough observations, and consequently the dashed arrows represent the constraints on the 
parameters of a and jY 1 . 



Obviously, for a — > oo and thus rj —> 1 which is the case considered in the WMAP analyses, 
the value of Zf approaches to its upper limit 10.7. In contrast, for a more realistic value of a 
on the order of unity, eq. (3.4) clearly shows that the redshift for full reionization could be 
much lower than its upper limit. 



4 Discussions 



The results presented above demonstrate that the effective power-law approximation for the 
high-redshift star formation history can dramatically simplify the calculations of reionization 
and CMB optical depth. More usefully, the connections among the star formation, reioniza- 
tion, and CMB optical depth could be approximately analytically addressed. These results 
make it convenient to test the theoretical models of star formation and reionization quickly 
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Figure 3. Numerically calculated evolutions of the ionization degree of hydrogen with redshift for 
different a (solid lines), where the values of jV 1 are taken according to eq. (3.1) with x(zh) = 0. The 
dashed lines represent power-law fittings to the numerical results. 

and effectively. At the same time, our results also provide an observational base to some 
relevant astrophysical studies, e.g., the estimation of the event rate of high-redshift GRBs. 

As implied by the probes to the Gunn-Peterson trough (Lya absorption) toward high- 
redshift quasars and galaxies, the full reionization of hydrogens could take place not far 
beyond zqp ~ 6 [37-40] . Of course, by considering that the Gunn-Peterson test is actually 
sensitive to small neutral fractions, the value of Zf could be somewhat higher. For example, 
detections of Lya-emitting galaxies at z = 6.6 [41, 42] suggested that the redshift for full 
reionization might be as high as z ~ 7. Therefore, for an independent and supplementary 
consideration, if we take 6 < Zf < 7, a constraint on the star formation history as well as on 
the parameter jKy could be derived from the JY^ — a and a — zj relationships: 

0.3 ± 0.2 < a < 1.3 ± 0.3, (4.1) 

and 

80 ± 12 < JCf < 230 ± 60. (4.2) 

On one hand, the range of a is well consistent with the one constrained by the observed near 
infrared excess in the cosmic infrared background spectrum [26]. On the other hand, for a 
typical value of 4000 of iV 7 for normal stars, the range of .JV^ indicates that the escaping 
fraction / esc is within the range from 2.0 ± 0.3% to 5.8 ± 1.5%. This is well consistent with 
the previous observations and calculations. 

The possible star formation history with 0.3 < a < 1.3 (shaded band) is presented in 
figure 1 in comparison with some observations. Then the comparison shows that, in general, 
the SFRs inferred from both GRBs and LBGs are too low to reionize the universe. To 
be specific, the slight shortage of the GRB-inferred SFRs may suggest that the calibration 
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Figure 4. The probability distributions of some cosmological parameters resulting from a fitting to 
the CMB power spectra with a hypcrbolic-tangent-function rcionization (dashed lines) or a power-law 
rcionization as eq. (3.2) (solid lines). 

between the GRBs and stars needs to be investigated more carefully. For example, some 
selection effects due to redshift measurements [43] should be taken into account, which are 
missed by ref. [4]. On the other hand, as also argued by some previous works (e.g., [4, 44, 45]), 
the significant discrepancy between the LBG-inferred SFRs and the required levels indicates 
that most of the high-redshift star formation could take place in very small halos hosting 
dwarf galaxies, which are too faint to be observed in deep LBG surveys. 

The result present above robustly exhibits the importance and advantage of the GRB 
observations in the SFR determination. In some previous works (e.g., [46]), the SFRs inferred 
from high-redshift GRBs were considered to be too high in contrast to the relatively low LBG- 
inferred SFRs. Therefore, Cheng et al. [46] proposed that these GRB-inferred SFRs may 
be overestimated due to that some high-redshift GRBs do not originate from the collapse 
of massive stars, but could be produced by the bursts of superconducting cosmic strings. 
However, the result addressed in this paper indicates an opposite situation, i.e., that the 
GRB-inferred SFRs actually are not overestimated but somewhat underestimated. Therefore, 
the expected number of the bursts of superconducting cosmic strings could be much lower 
than the one estimated in ref. [46] . 

Finally, it is worth noticing that the values of the cosmological parameters we adopted 
above actually are obtained by fitting the CMB power spectra with a hyperbolic-tangent- 
function reionization. Therefore, as pointed out by some previous works (e.g., [47, 48]), it 
seems necessary to revisit the fitting to the CMB power spectra with a more realistic reion- 
ization history, e.g., the power-law behavior as presented in eq. (3.2). For a preliminary 
attempt, with the CosmoMC code [49] (updated CAMB), we plot the probability distribu- 
tions of several cosmological parameters in figure 4, both for the cases of tangent-function 
and power-law reionizations. In comparison, the modification of the parameters due to the 
new reionization module is basically apparent and, specifically, the peak values of the pa- 
rameters shift about few percent. Nevertheless, such shifts would not significantly influence 
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the results obtained in this paper. 
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